Consistent initial data for CMD perturbations 
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Abstract 

We investigate the initial condition which are simultaneously consistent with per- 
turbation equations in both the radiation-dominated and the mater-dominated epochs. 
The exact formula for the spectrum transfer-function is derived. 

1 Introduction 

The scale invariance of the perturbations growth in the matter dominated epoch is known as 
the result of both the Newtonian and Einstein gravity theories. It is caused by the absence of 
the Euler term in the propagation equations [HE] - Partial differential equations derived as the 
first order expansion of Einstein equations, immediately reduce to ordinary ones. According 
to standard procedures the initial conditions for these equations are being imposed at the 
beginning of the matter dominated epoch. Initial velocity field is set to zero throughout 
the space, or alternatively, the decaying mode is neglected [3]. The perturbation then is 
the product of time-dependent factor and an arbitrary function of the space coordinates. 
Perturbations of any shape and any scale grow equally slow [H IU E] - too slow to form 
compact objects in the time interval admitted by observations [6j [7J |8]. 

Similar constraints for initial conditions are often applied to the n-body simulations [§]. 
Following [TD] numerous treatments additionally adopt the Zeldovich approximation [IT] or 
its generalizations p2]. The perturbations statistics and amplitudes are confronted with the 
CMB data and inflationary scenarios [131 E] • Most approaches adopt Newtonian gravity or 
the Newtonian limit of GR. The n-body simulations are promissory in non-linear regime. 
The question how to reach nonlinear regime still remains difficult for both simulations and 
analytical solutions. 
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In this paper we investigate initial condition which are simultaneously consistent with the 
perturbation equations in both, the radiation- and mater-dominated epochs. Following the 
techniques elaborated in the gravitational waves theory [TSJ [TBI El HBJ EE [20], we treat both 
epochs (with the perturbation corrections included) as the single Riemannian manifold. Ap- 
propriate junction conditions guarantee that the first and the second fundamental forms at 
the decoupling hypersurface are continuos. The perturbation modes are unique and reg- 
ular in the entire spacetime. The momentum spectrum of sound in the radiation-dominated 
epoch (the gauge-invariant and time-independent Fourier coefficients Ak) uniquely determine 
the perturbations evolution, both prior and past to decoupling. Due to junction conditions 
the perturbations in the late epoch (p = 0) inherit some properties from the epochs before 
decoupling. One of the symptoms is the breakdown of the late-time scale- invariance. 



2 Brief description of the method 



Einstein equations and Darmois-Israel conditions provide solution for the background scale 
factor 
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rjd stands for the time of decoupling, Ai is a constant of motion, Ai = a e 4 in the radiation 
filled universe. In the conformal Newtonian gauge with the metric tensor 



gfMU = a 2 ( V ) diag (-1 - 20, 1 - 20, 1 - 20, 1 - 20) . 



the propagation equations for the density perturbations in both epochs take forms 
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for the matter-dominated one. 

New perturbation variable \& defined as the Darboux transform 
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of the potential 
(6) 



allows to reduce Q to d'Alembert equation^] 
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V(t),x) = 0. 
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\l/ is an analogue of the massless scalar field. The general solution to ([7]) is the linear com- 
bination of running waves 

frfa x) = ^ Aufefa aQ + cc, (8) 



Ufe (77, a;) = exp — wr/)] 



where the modes 

ui,(n,x) = — 

are harmonic in rj, and normalized 

(u fc (r/,x), u fc / (77,2c)) = w(fc) 5 3 (k - k') 
according to the Klein-Gordon scalar product 

(ttifa x), $2(77, x)) = i f\ [*xd^* 2 - d<J». 
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The integral in formula (11) runs over arbitrary space-like hypersurface a 11 , therefore, the 



coefficients Ak are gauge-invariant and constant in time. Now, the field in the radiation- 
dominated epoch is a combination 



0(i) (v, x ) = Y1 AkU *w (n> x ) + c - c - 
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of modes 
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exp [i(kx — oorj)} . 
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In the matter dominated universe equation ^ is solved by an arbitrary combination of the 
constant or decreasing modes 
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1 For other gauges or gauge-invariant variables the equivalent procedures are given in 
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To assure the continuity of the first and the second fundamental forms at the decoupling 
hypersurface the coefficients Ak, &k and (3^ must be related by 
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Then the generic solution after decoupling is 

0(9 = $Z A^ (2) (rj, x) + c.c. 



where 
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The pair of solutions (13) and (18) combines into a single mode , w^ (2) } regular every- 
where in the space-time. The density contrast and the four-velocity corrections read 
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with <p given by ( 12 )— ( 13 ) or ( 17)— ( 18 ) respectively. The freedom to choose Ak is the only 
freedom left. Initial data on determined by ( 15 )— ( 16 ) are consistent with both equations, 



(Q and ([5]). While evolved forward or backward in time they restore the perturbation field 
in each of considered epochs. 



3 The transfer function 

We introduce dimensionless variables k, % and r based on the characteristic time-scale rjd 
(time of decoupling) and a characteristic length-scale rjd/y/3 (length of the sound horizon) 

K =^=k, X=^^, r = ^. (21) 
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k and r relate to observables 
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where H stands for the Hubble constant, z d is the redshift of decoupling, M. is the mass of 
sphere of the wavelength diameter. 

The density perturbations expressed in dimensionless variables take the form 



= ^2 AkUs m (r, x) + c - c - 

k 

6 m = Yl A kU 5(2) (r, x) + cc. 
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where 
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The linear dispersion relation co = fc/v3 is fulfilled. Perturbations of both the energy density 
and the expansion rate are continuous at the transition. 

The transfer function [2H] normalized to the decoupling epoch r = 1 reads 
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On the strength of (|27|) one has 
T 2 (k,t) 
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and therefore, 

a fc (r)4(r) = T 2 (k, r)a k (l)a* k (l) 
with the usual definition of a k given by 

k 

= ^ Qfc(r) exp(iKx) + c.c. 



(30) 
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The transfer function (29) differentiate perturbations of various length scales. The contour 



plot of T(Ai,z) as the function of mass and redshift (in the logarithmic scale) is shown 
on Figure 1. The mass M. is expressed in the solar mass units, the decoupling redshift 
Zd = 1200 is assumed. Conditions (15) and (16) do not restrain the perturbations spectrum. 
The invariant spectrum A k Al (or equivalently the coefficients ak{td) on S^) remains a free 
parameter of the model. 
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Figure 1: Logarithmic contour map of the transfer function lnT(lnA^,lnz). z stands for 
redshift, M. is the perturbation mass expressed in the Solar mass units. 
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